Cluster radioactivity in very heavy nuclei: a new perspective 
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Exotic cluster decay of very heavy nuclei is studied using the microscopic nuclear potentials 
obtained by folding density dependent M3Y effective interaction with the densities of the cluster 
and the daughter nuclei. The microscopic nuclear potential, Coulomb interaction and the centrifugal 
barrier arising out of spin-parity conservation are used to obtain the potential between the cluster 
and the daughter nuclei. Half life values are calculated in the WKB framework and the preformation 
factors are extracted. The latter values are seen to have only a very weak dependence on the mass 
of the emitted cluster. 
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The exotic phenomena cluster radioactivity in very 
heavy nuclei was suggested by Sandulescu et al. [l[ in 
1980. Since the first experimental observation of cluster 
radioactivity 0] , a lot of efforts have gone into the under- 
standing of the physics of cluster radioactivity. Predic- 
tions for a and various exotic decays have been made by 
the analytical superasymmetric fission model (ASAFM) 
[H, H[ with reasonable success. This was followed by the 
preformed cluster model (PCM) calculations which is dis- 
tinguished by the inclusion of the cluster preformation 
probability and was applied to a decay 0] with similar 
success. But both the theoretical approaches described 
above use phcnomcnological potentials for the nucleus- 
nucleus interactions. The ASAFM uses a parabolic po- 
tential approximation for the nuclear interaction within 
a superasymmetric fission model description which yields 
analytical expressions for the decay lifetimes, while the 
PCM uses a cos-hyperbolic form for nuclear interaction 
potential arising from the folding integral with suffi- 
ciently well known charge densities and the resulting nu- 
clear potential could be more or less accurately parame- 
terized to such a form. 

In the present work the phenomena of nuclear cluster 
radioactivity is studied theoretically using microscopic 
potentials within WKB framework of quantum tunnel- 
ing. The microscopic nuclear interaction potential is cal- 
culated by folding the density distributions of the emit- 
ted and daughter nuclei with density dependent M3Y 
effective interaction (DDM3Y). The DDM3Y effective 
nucleon-nucleon (NN) interaction used here was pevi- 
ously used successfully for elastic and inelstic scattering 
of protons M, 0] , proton radioactivity [|[ and a radioac- 
tivity [§, [l(| whose density dependence was obtained 
from nuclear matter calculations [ll|. The cluster pre- 
formation factors are extracted from the calculated and 



the measured half lives of cluster radoactivity and its 
systematics are studied. 

The decay constant A for cluster radioactivity is a 
product of cluster preformation probability Pq in the 
ground state, the tunneling probability through barrier 
P and the assault frequency v. Cluster emission half 
life T1/2 = In 2/ A calculated within WKB framework for 
barrier penetration probability P is: 



T 1/2 = [(hln2)/(2E v )][l+exp(K)} 



(1) 



where E v , the zero point vibration energy, is assumed 
to be proportional to Q-value of the spontaneous cluster 
emission. The WKB action integral is given by 



K = 



y/2n(E(R) - E v - Q)dR 



(2) 
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where the total interaction energy between the emitted 
cluster and the daughter nucleus E(R)=Vn{R)+Vc(R) + 
h 2 l(l + l)/(2fiR 2 ) is equal to the sum of the nuclear in- 
teraction energy, Coulomb interaction energy and the 
centrifugal barrier arising from spin-parity conservation. 
The reduced mass fj, = mA e Ad/ A where A e ,Ad,A are 
the mass numbers of the emitted cluster, the daughter 
and the parent nuclei respectively and m is the nucleonic 
mass measured in the units of MeV /c 2 . R a and Rb are the 
second and third turning points of the WKB action inte- 
gral determined from the condition: E(R a ) = Q + E v — 
E(Rb) which provides three classical turning points. The 
energetics allow spontaneous emission of clusters only if 
the released energy Q = M — (M e + M4) is a positive 
quantity, where M, M e and Md are the atomic masses (in 
MeV) of the parent nucleus, the emitted cluster and the 
residual daughter nucleus, respectively. The microscopic 
nuclear potentials Vjv(iZ) are obtained by double fold- 
ing the densities of the emitted cluster and the residual 
daughter nucleus with the finite range realistic density 
dependent M3Y effective interacion as 
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V N (R) =/ / Pi(rl)p 2 (r5)«( S )d 3 r 1 d 3 r 2 (3) 

where s=|r2 — rf + R\, p\ and p2 are the density dis- 
tributions for the two composite nuclear fragments hav- 
ing the spherically symmetric form given by p\p[f) = 
po/[l + exp((r — c)/a)\ where c = r p (l — n 2 a 2 /3r 2 ), r p 

1 /3 

= 1.13A d ' e , a — 0.54 fm and the value of p is fixed by 
equating the volume integral of the density distribution 
function to the mass number of the fragment. The exper- 
imental charge density distributions in case of the heavier 
nuclei can be well described by the two parameter Fermi 
function [l2| and since the charge i. e. the proton (p) and 
the neutron (n) density distributions should have similar 
forms due to the same strengths of the n-n and p-p nu- 
clear forces, the matter density distribution can be well 
described by the spherically symmetric Fermi function. 
The density distribution function in case of a particle 
has the Gaussian form p 2 (r) = 0.4229 exp(-0.7024r 2 ) 
whose volume integral equals A a (= 4), the mass number 
of the a particle. 

The density dependent M3Y [l3| effective interaction 
(DDM3Y) [14| appearing in Eq.(3) is given by v(s) — 
t M3Y (s,E)g( Pl ,p 2 ) where 

t Msr = 7999 ex P(~ 4s ) _ 21M ^~^1 
(4s) (2.5s) 

-276(1 - 0.005E/A e )S{s) (4) 

where E and A e are the laboratory energy and projec- 
tile mass number respectively. In the present case of 
radioactivity, it can be shown that E/A e = Qm/ pi where 
m and p are the nucleonic mass and reduced mass of 
the A e + Ad system, respectively, in units of MeV/c 2 . 
The zero-range potential represents the single-nucleon 
exchange term and the density dependence g(pi,P2) — 
C(l — /3pi^ 3 )(l — (3p^ 3 ) takes care of the higher order ex- 
change effects and the Pauli blocking effects. This density 
dependent M3Y effective NN interaction supplemented 
by the zero-range potential is used to determine the nu- 
clear matter equation of state. The equilibrium density 
of the nuclear matter is determined by minimizing the 
energy per nucleon. The density dependence parame- 
ters have been fixed by reproducing the saturation en- 
ergy per nucleon and the saturation density of spin and 
isospin symmetric cold infinite nuclear matter. Although 
the density dependence parameters for single folding can 
be determined from the nuclear matter calculations and 
used successfully for proton radioactivity [1] and scatter- 
ing 0, 0] , the transition to double folding is not straight- 
forward. The parameter /3 can be related to mean free 
path in nuclear medium, hence its value should remain 
same ~ 1.6/m 2 as obtained from nuclear matter calcula- 
tions [llj while the other constant C which is basically an 
overall normalisation constant may change. The value of 
this overall normalisation constant has been kept equal 



to unity which has been found ~ 1 from optimum fit to 
a large number of alpha decay lifetimes. The density- 
dependence of the effective projectile-nucleon interaction 
is found to be fairly independent of the projectile, as 
long as the projectile-nucleus interaction is amenable to 
a single-folding prescription. This argument can be fur- 
ther stretched to mean that, in a double folding model, 
the density-dependent effects on the nucleon-nucleon in- 
teraction can be factorized into a 'target term' times a 
'projectile term' (l5j . 

The cluster preformation probability Pq for a is close to 
unity, e.g. 1, 0.6, 0.35 for even-even, odd-even and odd- 
odd emitters @ respectively. The results of the present 
theoretical calculations for few a emitters are listed in 
Table I. Same set of a and cluster emitters of reference 
fill are chosen. The zero point vibration energies used in 
the present calculations are the same as that described in 
reference pj} immediately after eqn.(4) for the a cluster 
and by eqns.(5) for the heavier clusters. The shell effects 
for every cluster radioactivity are implicitly contained in 
the zero point vibration energy due to its proportionality 
with the Q value, which is maximum when the daugh- 
ter nucleus has a magic number of neutrons and protons. 
Values of the proportionality constants of E v with Q is 
the largest for even-even parent and the smallest for the 
odd-odd one. Other conditions remaining same one may 
observe that with greater value of E v , lifetime is short- 
ened indicating higher emission rate. 



TABLE I: Half lives of a decay obtained in the present cal- 
culation. The asterisk symbol represents theoretical values 
assuming a preformation probability to be unity. 



Parent 


Cluster 


Q (MeV) 


logT ex 


logTth 


Po 


logT th 


/ 


212 Po 


4 He 


8.954 


-6.52 


-7.01 


1.00 


-7.01 





213p 


4 He 


8.536 


-5.44 


-5.65 


0.60 


-5.43 





214 Po 


4 He 


7.833 


-3.78 


-4.06 


1.00 


-4.06 





215 At 


4 He 


8.178 


-4.00 


-4.47 


0.35 


-4.01 






The theoretical half lives of cluster radioactivity for 
very heavy nuclei are calculated assuming cluster pre- 
formation factor to be unity. Hence the preformation 
factors Po can be calculated as the ratios of the calcu- 
lated half lives to the experimentally observed half lives 
and listed in Table II. This represents the preformation 
factor and may be considered as the overlap of the actual 
ground state configuration and the configuration repre- 
senting the cluster coupled to the ground state of the 
daughter. Po for normal a emitters is close to unity 
Q . Superheavy emitters being loosely bound than highly 
bound a, Pq is expected to be high and present calcula- 
tions with Pq=1 provide excellent description [jj, [ToL Il8j 
of a decay for recently discovered superheavy nuclei. For 
weakly bound heavy cluster decay it is expected to be 
orders of magnitude less than unity. 

The basic problem of the present theoretical study is 
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TABLE II: Half lives assuming preformation factor to be unity 
and corresponding preformation factors of cluster decay ob- 
tained in the present calculation. The asterisk symbol rep- 
resents theoretical values assuming preformation factor to be 
unity and using normalisation of 0.7 for the nuclear potentials. 

Parent Cluster Q (MeV) logT ex logT* h logTth -logPo I 

221 Fr 14 C 31.317 14.52 13.46 11.71 2.81 3 

221 Ra 14 C 32.396 13.39 13.23 11.47 1.92 3 

222 Ra 14 C 33.050 11.00 10.38 8.71 2.29 

223 Ra 14 C 31.829 15.20 14.34 12.56 2.64 4 

224 Ra 14 C 30.540 15.92 15.12 13.34 2.58 

225 Ac 14 C 30.477 17.34 17.25 15.39 1.95 4 

226 Ra 14 C 28.200 21.34 20.22 18.35 2.99 

228 Th 20q 44.720 20.72 21.02 18.90 1.82 

230 U 22 Ne 61.400 19.57 20.81 18.23 1.34 

230 Th 24 Ne 57.57! 24.64 25.07 22.58 2.06 

231 Pa 24 Ne 60.417 23.38 22.74 20.32 3.06 1 

232 U 24 Ne 62.310 20.40 20.67 18.30 2.10 

233 U 24 Ne 60.486 24.82 24.31 21.82 3.00 2 

234 U 24 Ne 58.826 25.25 25.97 23.38 1.87 

233 U 25 Ne 60.776 24.82 24.35 21.92 2.90 2 

234 U 26 Ne 59.466 25.07 26.27 23.84 1.23 

234 U 28 Mg 74.110 25.74 25.66 22.93 2.81 

236p u 28 Mg 79.670 21.67 21.38 18.80 2.87 

238p u 28 Mg 75.912 25.70 26.25 23.49 2.21 

238 Pu so Mg 76.824 25.28 25.85 23.22 2.06 

238p u 32 gi 9i.i9o 25.30 26.21 23.30 2.00 

242 Cm 34 Si 96.509 23.15 23.82 21.16 1.99 



that there is no guarantee that the experimentally ob- 
served decays proceed from the ground state of the parent 
nucleus to that of the daughter nucleus which is assumed 
in the present calculations. It is a fundamental difficulty 
associated with the decay of the odd mass parent nu- 
clei or the odd mass emitted clusters. These nuclei may 
be decaying predominantly to a low-lying excited state of 
the daughter nucleus. When the exotic cluster is removed 
from the parent nucleus, the state of the core left over 
may be quite different from that of the ground state of the 
daughter nucleus, but rather similar to that of one of its 
excited states. Hence, decays would go preferentially to 
this particular excited state, and be strongly suppressed 
to the ground (and other) states. Present model is not 
microscopic enough to predict, a priori, which excited 
state may be most appropriate for this role. However, 
since the effects due to nuclear structure do not appear 
explicitly, present calculations can also provide the the- 
oretical half lives for transitions other than ground state 
to ground state if only the spin-parities and the corre- 
sponding experimental Q values are precisely known. 

A systematic study of even-even parent and daughter 
combinations for which centrifugal barrier is zero, it was 
shown [19| that —logPg can be fitted either to c\ZdZ e +C2 



or to c-iA e +Ci where Zd, Z e are the charge numbers of the 
daughter and emitted nuclei. The folded nucleus- nucleus 
interaction potentials were obtained using M3Y effective 
interaction without any density dependence and the re- 
sulting nuclear potentials were renormalized (RM3Y) by 
multiplying them with a factor of 0.55 [l||. It is worth- 
while to mention here that due to attractive character 
of the M3Y forces the saturation conditions for cold nu- 
clear matter is not fulfilled [20]. The realistic descrip- 
tion of nuclear matter properties can be obtained only 
with the density dependent M3Y effective interaction. 
Using DDM3Y effective interaction the saturation condi- 
tions can be fulfilled and in turn the constants of density 
dependence can be obtained from the saturation condi- 
tions of the cold symmetric nuclear matter. Hence, in 
the present case the nuclear potentials are not renor- 
malised but the same density dependent effective interac- 
tion DDM3Y is used whose density dependence was ob- 
tained from nuclear matter calculations based on Hartree 
or mean field assumption [ll[ and was used quite suc- 
cessfully for proton scattering [a, [7J , proton radioactivity 
Q and a radioactivity [j| [ToL Il8|^ Present calculations 
yield constants c\ = -0.38 x 10 -3 , C2 = 2.56 from a very 
poor quality of fit and c 3 = 0.028, c 4 = 2.10 from some- 
what better fit, which are quite different from reference 
These results imply absolutely no dependence on 
ZdZ e and very weak dependence on A e . The preforma- 
tion probabilities found for the range of clusters consid- 
ered here are « 10~ 2 to 10 -3 which tend to disagree with 
the results of a potential model with preformed clusters 
[Hi HH H3 ■ It is interesting to see in Table I that for a 
cluster the agreement is excellent. Although present cal- 
culations provide excllcnt estimates for the a decay half 
lives without adjusting the depth of the nuclear poten- 
tials, it is interesting to note that the theoretical cluster 
radioactivity half lives {T^ h ) are in reasonable agreement 
with experimental data spanning about fifteen orders of 
magnitude if a renormalisation of 0.7 for the nuclear po- 
tentials is used. 

To summarise, the half lives for cluster radioactivity 
have been calculated with microscopic nuclear potentials 
which are based on profound theoretical basis. Present 
calculation is found to be highly successful for a radioac- 
tivity including that of superheavies. For cluster radioac- 
tivity, present calculations with preformation probability 
Pq = 1 systematically underestimate half lives, typically 
by factors of 10~ 2 to 10 -3 . The cluster preformation 
factors Po extracted from the ratio of calculated to mea- 
sured half lives of cluster radoactivity shows very weak 
dependence on the mass number of the emitted nucleus 
and is roughly constant over the limited range of clusters 
(6< Z e <14) considered here. The implication of such 
a state of affairs suggests contrary to all previous expec- 
tations, that it may be possible to extract nuclear struc- 
ture information from exotic decay data. There seems 
to be a distinct chance that the single-particle orbitals 
filled by the daughter (core) nucleons are distorted and 
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change their ordering due to polarisation effects caused 
by the presence of a large cluster. This, in turn, in- 
fluences which excited state of the daughter nucleus is 
predominantly populated. It is interesting to note that 
the theoretical cluster radioactivity half lives (T^ h ) are in 
reasonable agreement with experimental data spanning 
about fifteen orders of magnitude if a renormalisation of 
0.7 for the microscopic nuclear potentials is used. Many 



questions remain to be answered concerning the univer- 
sality of this phenomenon, and whether different clusters 
cause varying amounts of core polarisation, so that dif- 
ferent excited states of the same daughter nucleus are 
preferred when different clusters are emitted, and offers 
a very interesting challenge to more microscopic models, 
and to the experimentalists. 
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